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1. INTRODUCTION 

Today’s power system, comprising of generator, transmission lines and different types of loads is 
a complex network. As the power demand increase some transmission line is more loaded than was planned. 
In order to meet the growing power demand, utilities have an interest in better utilization of available power 
system capacities, existing generation and existing power transmission network, instead of building new 
transmission lines and expanding substations. The main drawback of the power system is deteriorating 
voltage profiles and decreasing system stability and security; the main cause of this drawbacks are due to 
the overloaded of the power flow in transmission lines of power system. In addition, existing ancient 
transmission facilities, in most cases, are not designed to handle the control requirements of complex and 
highly interconnected power systems. This overall scenario force us to review of ancient transmission 
methods and practices, and the creation of new concepts, which might allow the use of existing generation 
and transmission lines up to their full capabilities without reduction in system stability and security [1]. 
The main parameters which are responsible to determine the transmitted electrical power over a line are line 
impedance, the receiving and sending ends voltages, and phase angle between the voltages. Therefore, by 
controlling, one or more of the transmitted power factors; it is possible to control the active as well as 
the reactive power flow over a line. Small-signal stability is the ability of the system to return in a normal. 

Operating state following a small disturbance. Investigations involving this stability concept usually 
involve the analysis of the linearized state space equations that define the power system dynamics. To damp 
out the low power oscillation frequency and increase system oscillations stability, the installation of Power 
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System Stabilizer (PSS) is both economical and effective [2]. Recently appeared FACTS (Flexible AC 
Transmission System)-based stabilizer offer an alternative way in damping power system oscillation. 
The primary function of the FACTS controllers is not only it Damping Duty, but also to increase the overall 
power system oscillation damping characteristics [3]. 

Many researchers have presented dynamic models of UPFC [4, 5-14]. A modified linearised 
Phillips-Heffron model of a power system installed with UPFC have presented by Wang [4, 13-14]. 
The objective of this paper is to design a UPFC based Power Oscillation Damping (POD) controller to damp 
the low frequency electromechanical oscillations over wide range of operating conditions [4]. The procedure 
to achieve the paper objective is as follows: 

— The model (Phillips—Heffron) of the single machine infinite bus (SMIB) power system installed with 
UPFC is obtained by linearising the non-linear equations around a nominal operating point [4]. 

— To present systematic approach for designing UPFC based power oscillation damping controller. 

— Design a POD controller using LQR technique which places the eigenvalue corresponding to mode of 
oscillation at desired location such that eigenvalues get placed within a vertical degree of stability. 

— To demonstrate the effectiveness of the designed POD controller under different controlling parameter. 


2. INVESTIGATED SYSTEM 

Figure 1 shows a single machine infinite bus power system installed with unified power flow 
controller between bus A and bus B on the transmission line, which consists of an excitation transformer 
(ET), a boosting transformer (BT), two three-phase GTO based voltage source converters (VSC’s), a de link 
capacitor and based on pulse width modulation converters (assumed). The system parameters are given in 
Appendix-A. 





Figure 1. UPFC installed in SMIB power system 


3. UNIFIED POWER FLOW CONTROLLER 

The dc link provides a bi-directional path to exchange active power between the converters, coupled 
a static synchronous compensator and static synchronous series compensator. me, mb , de, db are the input 
control parameter of the UPFC. Here me is amplitude modulation ratio of shunt VSC-E, de is phase angle of 
shunt VSC-E, mb is amplitude modulation ratio of series VSC-B and ôb represents phase angle of series 
VSC-B [15, 16]. Power system which comprises a synchronous generator connected to an infinite bus 
through a transmission line and stepping up transformer. The generator is assumed to have Automatic 
Voltage Regulator (AVR) controlling its terminal voltage. The UPFC is used in this study to just analyze 
the purpose of power system stability and its characteristics. The UPFC can fulfill multiple control 
objectives. 


4. MODIFIED PHILLIPS-HEFFRON MODEL OF A SMIB SYSTEM WITH UPFC 

Figure 2 shows the Modified Phillips-Heffron model of power system installed with UPFC, developed 
by Wang [4] with the modification of the basic Phillips-Heffron model including UPFC. Around a nominal 
operating point, by linearising the nonlinear model, this model has been developed. 

The parameters of the model depend on the system parameters and the operating condition. In this 
model [Au] is the column vector while [Wpu], [Wau], [Wvu] and [Weu] are the row vectors. 
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Where, 


[Au] = [Am Ade Am Ady] 





[Wu] 
Wal = 
[Wv] = 


[Weu] = 


= [Woe Woe Wop Wa] 
[We Wade Wo Waa] 
[Wye Wyae Ww Wvo] 


[Wee Wede Wo Wea] 


The control parameters of the UPFC are, 
m»—pulse width modulation index of series inverter. The magnitude of series injected voltage can be 
controlled by controlling mp. 

—  6b-—Phase angle of series inverter. 6b when controlled results in the real power exchange. 

— me-—pulse width modulation index of shunt inverter. By controlling me, the voltage at a bus where UPFC 
is installed, is controlled through reactive power compensation. 


—  6e-Phase angle of the shunt inverter, which regulates the dc voltage at dc link. 
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Figure 2. Modified Phillips-Heffron model of power system installed with UPFC 


5. MODEL ANALYSIS 
The constant parameters of the model computed for nominal operating condition and system 


parameters are. 


W1 = 0.1372 
W2 = 0.4350 
= 0.4727 
W4 = 0.0598 
Ws = -0.0159 
We = 0.5092 
W7 = 80.9318 
Ws = 21.0677 
Wo = 40.6634 


Int J Rob & Autom, Vol. 


Wpb =0.0715  Wpae = -0.2514 
Wap = -0.0223 Waa = - 0.2243 
Ww = 0.0145 Wve = 0.0782 
Wye = 0.7860 Wen = 0.1763 
Wee = -0.2451 Wee = 0.0018 
Wye =0.1597 Wew = 37.7306 
Wpa = -0.0229  Wede = 77.2952 
Waa = - 0.0204 Wpa = 0.4287 
Wæ = 0.0061 Waa = -0.1337 
Wya= 0.0871 
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6. DESIGN OF POD CONTROLLER (LINEAR QUADRATIC REGULATOR TECHNIQUE) 
In Figure 3 is generalized Diagram of POD Controller. The linearized state-space model of SMIB 
power system is obtained by phillip-heffrons model as expressed by (1), 


X = AX + BU 0) 


Where A and B are the matrices of the system and input respectively. X is the system state vector, 
and U is the input state vector. The matrices A and B are constant under the assumption of system linearity. 
Determine the matrix K of the LQR vector, that is, if we set u(t)= -Kx(t), so in order to determine 
the performance index (2) [17], [18]. 


J = | (x7 Qx +u" Ru) dt (2) 


In as shown in (48) the second term on the right side accounts for the expenditure of the energy on 
the control efforts. The matrix Q and R determine the relative importance of the error and the expenditure of 
this energy. The as shown in (48) can be rewrite as (3) 


œ 


J= J (x Ox+ x" K"RKx) dt =| x’ (Q+ K" KR)xdt (3) 


10) 


For solving the parameter-optimization problem, set the equation as shown in (4-5) 


x’ (O+K'RK)x= -Z (x" Px) (4) 
t 
x’ (QO+K'RK)x =—x' Px—x" Px (5) 
=—x"[(A— BK)’ P+ P(A-BK)]x 
Comparing both sides of the as shown in (5) and this equation must hold true for any x, (6) 
(A— BK)’ P+ P(A— BK) =—(O+ K’ RK) (6) 


Since R is a positive-definite Hermitian or real symmetric matrix (assume), then write it, (R=T'T) 
where T is a nonsingular matrix, and (7) 


A’ P+ PA+[TK —(T7 )'B’ P|— PBR"'B'’P+O=0 (7) 
The minimization of J with respect to K requires the minimization of (8) 

x [TK — (T7) B" PY [TK —(77 )' B P\x (8) 
As shown in (8) is nonnegative, the minimum occurs when it is zero, or when (9-10) 

TK =(T')'B'™P (9) 

K=T'(')'B’P=R'B'P (10) 
Thus a control law is (11) 

u(t) =—Kx(t) =—R™'B" Px(t) (11) 
In which P must satisfy the reduced Riccati equation (12) 


ATP+PA-PBR'BTP+Q=0 (12) 
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Steps for controller design 

— Required data—A, B, Q, R and N matrix 

— Data for calculation-K matrix, P matrix and eigenvalue “e” 

— Set N matrix is zero 

— QandR are the positive definite real symmetric matrix. 

— Calculate the P matrix with the help of the Riccati equation (12) 
— Find the value of K matrix from the equation (10) 
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Figure 3. Generalized diagram of POD controller 


7. SIMULATION RESULTS UNDER DIFFERENT SYSTEMS AND AT VARIOUS LOADING 
CONDITIONS 

The proposed model of UPFC in single machine infinite bus power system Figure 1 has been used in 
order to study the damping performance. To study the performance of the proposed controller, simulation 
results under different system conditions and at various loading conditions i.e. at normal operating point 
corresponding to line loading of 1.0 pu and at 20 percent decrease and increase in line loading are shown. 
It can be readily seen that the proposed controller performs better in terms of reduction of overshoot and 
settling time than system without UPFC and system with UPFC only. Simulation results with variation in 
system- state, rotor angle (5) of generator is only considered. 

In Figure 4 is variation in delta with UPFC, without UPFC and with POD controller of unstable 
(D=0) system at nominal loading (1.0 pu). In Figure 5 is variation in delta with UPFC, without UPFC and 
with POD controller of unstable (D=0) system at 20 percent decrease in loading (0.8 pu). In Figure 6 is 
variation in delta with UPFC, without UPFC and with POD controller of unstable (D=0) system at 20percent 
increase line loading (1.2 pu). 

In Figure 7 is variation in delta without UPFC, with UPFC and with POD controller of weak (D=4) 
system at nominal line loading (1.0 pu). In Figure 8 is Variation in delta with UPFC, without UPFC and with 
POD controller of weak (D=4) system at 20 percent increase line loading (1.2 pu). In Figure 9 is variation in 
delta with UPFC, without UPFC and with POD controller of stronger (D=8) system at nominal loading 
(1.0 pu). In Figure 10 is variation in delta with UPFC, without UPFC and with POD controller of strong 
(D=8) system at 20 percent decrease line loading (0.8 pu). 
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Figure 4. Variation in delta of unstable (D=0) system at nominal loading (1.0 pu) 
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In Figure 11 is variation in delta with UPFC, without UPFC and with POD controller of strong 


(D=8) system at 20 percent increase line loading (1.2 pu). In Figure 12 is variation in delta with UPFC, 
without UPFC and with POD controller of weak (D=4) system without any line loading. In Figure 13 is 
variation in delta with UPFC, without UPFC and with POD controller for modulation index of both 
converters 0.6. In Figure 14 is variation in delta with UPFC, without UPFC and with POD controller for both 


converters phase angle 0.6. 
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Figure 5. Variation in delta of unstable (D=0) system at 20 percent decrease in loading (0.8 pu) 


Variation in delta of SMIB system w ithout UPFC,w ith UPFC and w ithcontroller UPFC 
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Figure 6. Variation in delta of unstable (D=0) system at 20percent increase line loading (1.2 pu) 
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Figure 7. Variation in delta of weak (D=4) system at nominal line loading (1.0 pu) 
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Variation in delta of SMIB system w ithout UPFC,w ith UPFC and w ithcontroller UPFC 
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Figure 8. Variation in delta of weak (D=4) system at 20 percent increase line loading (1.2 pu) 
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Figure 9. Variation in delta of stronger (D=8) system at nominal loading (1.0 pu) 


Variation in delta of SMIB system w ithout UPFC,w ith UPFC and w ithcontroller UPFC 





























4 r r r r r r r r r 
with UPFC controller with UPFC 
\ without UPFC 
Z | 4 
5 Y | 
= TE 
aJ t | | | \ = 
eal YY d 
“3h IR 4 
-4 r r r r r r r r r 
o 1 2 3 4 5 6 7 8 9 10 


Time(sec) (sec) 


Figure 10. Variation in delta of strong (D=8) system at 20 percent decrease line loading (0.8 pu) 


Int J Rob & Autom, Vol. 8, No. 3, September 2019 : 164 — 173 


Int J Rob & Autom ISSN: 2089-4856 o 171 





Variation in delta of SMIB system w ithout UPFC,w ith UPFC and w ithcontroller UPFC 
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Figure 11. Variation in delta of strong (D=8) system at 20 percent increase line loading (1.2 pu) 


Variation in delta of SMIB system w ithout UPFC,w ith UPFC and w ithcontroller UPFC 
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Figure 12. Variation in delta of weak (D=4) system without any line loading 
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Figure 13. Variation in delta for modulation index of both converters 0.6 
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Figure 14. Variation in delta for both converters phase angle 0.6 
EIGEN VALUE ANALYSIS 
Eigen values with UPFC POD controller for weak stable SMIB system as shown in Table 1. 
Table 1. Eigen values with UPFC POD controller for weak stable SMIB system 
(PF=.85, ME =.5, AE =.5, MB =.5, AB =.5) 
Controller f System with System with UPFC 
Damping constant (D) ee mie UPFC only controller 
-98.7309 -102.18 -102.18 
0.0224 + 5.61801 -33.02 -33.02 
D=0 0.0224 - 5.6180i -1.26 + 4.631 -3.79 + 4.631 
- 1.6663 -1.26 — 4.631 -3.79 — 4.631 
0 -3.35 -3.35 
-98.7309 -102.18 -102.18 
-0.2270 + 5.6116i -33.04 -33.04 
D=4 -0.2270 - 5.6116i -1.62 + 4.70i -4.86 + 4.70i 
-1.6674 -1.62 — 4.70i -4.86 — 4.70i 
0 -3.12 -3.12 
-98.7309 -102.18 -10218 
-0.4765+5939i -33.05 -33.05 
D=8 -0.4765-5.5939i -1.98 + 4.781 -5.95 + 4.781 
- 1.6686 -1.98 — 4.781 -5.95 — 4.78i 
0 -0.0289 -2.89 
CONCLUSION 


In this paper, the power system low frequency electromechanical oscillations was damped via LQR 
technique (using MATLAB tool) based POD controller when applied independently with UPFC and 
investigated for a SMIB power system. For the proposed controller design problem, an eigenvalue-based 
objective function to maximize the system damping ratio among all complex eigenvalues was developed. 
The effectiveness of the proposed controller in damping the low frequency EM mode of oscillations and 
hence improving power system dynamic stability have been verified. 
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